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Abstract 



In this paper, we derive Gallager's random coding error exponent for multiple-input multiple-output 
(MIMO) channels, assuming no channel-state information (CSI) at the transmitter and perfect CSI 
at the receiver. This measure gives insight into a fundamental tradeoff between the communication 
reliability and information rate of MIMO channels, enabling to determine the required codeword length 
to achieve a prescribed error probability at a given rate below the channel capacity. We quantify the 
effects of the number of antennas, channel coherence time, and spatial fading correlation on the MIMO 
exponent. In addition, general formulae for the ergodic capacity and the cutoff rate in the presence of 
spatial correlation are deduced from the exponent expressions. These formulae are applicable to arbitrary 
structures of transmit and receive correlation, encompassing all the previously known results as special 
cases of our expressions. 

Index Terms 

Block fading, channel capacity, cutoff rate, multiple-input multiple-output (MIMO) system, random 
coding error exponent, spatial fading correlation. 



The channel capacity is a crucial information-theoretic perspective that determines the funda- 
mental limit on achievable information rates over a communication channel [1]. However, since 
the channel capacity alone gives only the knowledge of the maximum achievable rate, a stronger 
form of the channel coding theorem has been pursued to determine the behavior of the error 
probability P c as a function of the codeword length iV and information rate R [2]-[4]. The 
reliability function or the error exponent of a communication system is defined by [2] 



where P° pt (R, N) is the average block error probability for the optimal block code of length N 
and rate R. 1 The error exponent describes a decaying rate in the error probability as a function 
of the codeword length, and hence serves to indicate how difficult it may be to achieve a certain 
level of reliability in communication at a rate below the channel capacity. Although it is difficult 
to find the exact error exponent, its classical lower bound is available due to Gallager [3]. This 

'in the following, we will use the term "error probability" to denote the average block error probability. 



I. Introduction 




N 



N^oo 
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lower bound is known as the random coding error exponent or Gallager's exponent in honor of 
his discovery, and has been used to estimate the codeword length required to achieve a prescribed 
error probability [5]-[7]. 

The random coding exponent was extensively studied for single-input single-output (SISO) and 
single-input multiple-output (SIMO) flat-fading channels with average or peak power constraint 
[5], [6]. For SIMO block-fading channels, the random coding exponent was derived in [8] with 
perfect channel-state information (CSI) at the receiver, where it has been shown that although 
the capacity is independent of the channel coherence time (first asserted in [9] and also recently 
addressed in [10] and [11] for multiple- antenna communication), the error exponent suffers a 
considerable decrease due to a reduction in the effective codeword length as the coherence 
time increases. 2 Therefore, this so-called channel-incurable effect reduces the communication 
reliability. While there are numerous prior investigations (following the seminal work of [12]- 
[15]) on the capacity for multiple-input multiple-output (MIMO) channels [16]— [23], only limited 
results are available for error exponents. The random coding exponent were given implicitly 
in [16] (without final analytical expressions) for independent and identically distributed (i.i.d.) 
Rayleigh-fading MIMO channels with a single-symbol coherence time, perfect receive CSI, and 
Gaussian inputs subject to the average power constraint. Also, the random coding exponent was 
analyzed in [7] for i.i.d. block-fading MIMO channels with no CSI and isotropically unitary 
inputs subject to the average power constraint. 

In this paper, taking into account spatial fading correlation, we derive Gallager's exponent 
for MIMO channels. We consider a block-fading channel with Gaussian inputs subject to the 
average power constraint and perfect CSI at the receiver. Our results resort to the methodology 
developed in [22] and [23], which is based on the finite random matrix theory [24], [25]. The 
MIMO exponent obtained in the paper provides insight into a fundamental tradeoff between 
the communication reliability and information rate (below the channel capacity), enabling to 
determine the required codeword length for a prescribed error probability. It is interesting 
to note that as a special case of this reliability-rate tradeoff, one can obtain the diversity- 
multiplexing tradeoff of MIMO channels [12], [13], [26], which is a scaled version of the 

2 This observation is parallel to the divergent behavior of the channel capacity and cutoff rate of a channel with block memory 
[9]. 
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asymptotic reliability-rate tradeoff at high signal-to-noise ratio (SNR). We quantify the effects of 
the number of antennas, the channel coherence time, and the amount of spatial fading correlation 
on the MIMO exponent. Moreover, the general formulae for the ergodic capacity and cutoff rate 
are deduced from the exponent expressions. In particular, our capacity formula embraces all 
the previously known results for i.i.d. [16], [21], one-sided correlated [19], [20], and doubly 
correlated [22] channels. 

This paper is organized as follows. In Section [ill signal and channel models are presented. 
SectionlrTTlderives the expression for the MIMO random coding exponent. Section |lV]gives proofs 
of the main results stated in Theorem d In Section |Vj some numerical results are provided 
to illustrate the reliability-rate tradeoff in block-fading MIMO channels. Finally, Section |VI| 
concludes the paper. 

Notation: Throughout the paper, we shall use the following notation. N, R, and C denote 
the natural numbers and the fields of real and complex numbers, respectively. The superscripts 
T and f stand for the transpose and transpose conjugate, respectively. I n is the n x n identity 
matrix and (Aij) denotes the matrix with the (i,j)th entry A^. The trace operator of a square 
matrix A is denoted by tr (A) and etr (A) = e tr( - A \ The Kronecker product of matrices is 
denoted by <g>. By A > 0, we denote A is positive definite. For a Hermitian matrix A E C nxn , 
Ai (A) > A 2 (A) > . . . > A n (A) denotes the eigenvalues of A in decreasing order and A (A) 6 
]R n denote the vector of the ordered eigenvalues, whose ith element is A^ (A). Also, g(A) 
denotes the number of distinct eigenvalues of A, and A^ (A) and \k (A), k = 1, 2, . . . , g (A), 
denote the distinct eigenvalues of A in decreasing order and its multiplicity, respectively, that is, 
A/^ (A) > X, 2 , (A) > . . . > A/ (A) and YltLi Xk (A) = n. Finally, we shall use the notation 
X E C mxn ~ M m , n (M. E. V) to denote that a random matrix X is (matrix-variate) Gaussian 
distributed with the probability density function (pdf) 

px (X) = u- mn det (E)" n det (*)""' etr { -E 1 [X - M) * 1 (X - M) 1 } (1) 

where M E C mxn , E = E f E C mxm > 0, and = & E C nxn > 0. 

II. Signal and Channel Models 

We consider a MIMO system with transmit and receive antennas, where the channel 
remains constant for N c symbol periods and changes independently to a new value for each 



December, 2005 



DRAFT 



4 



SUBMITTED TO THE IEEE TRANSACTIONS ON COMMUNICATIONS 



coherence time, i.e., every N c symbols. Since the propagation coefficients independently acquire 
new values for every coherence interval, the channel is memory less when considering a block 
length of N c symbols as one channel use with input and output signals of dimension ut x N c 
and n R x N c , respectively. 

For an observation interval of N b N c symbol periods, the received signal is a sequence {Y k }^ v 
each Y k E C" rxJVc is given by 

Y k = H k X k + W k , k = 1, 2, . . . , N h (2) 

where X k E C 1 ^ 1 are the transmitted signal matrices, H k E C nRXTlT are the channel matrices, 
and W k ~ M nn ,Nz (0, N Q I nn , I at c ) are the additive white Gaussian noise (AWGN) matrices. Fig.[T] 
shows a communication link with n T transmit and n R receive antennas to communicate at a rate 
R (in bits or nats per symbol) over N h independent iV c -symbol coherence intervals. Since the 
channel is memoryless with identical channel statistics for each coherence time interval, the 
index k can be dropped. 

Let px (X) be the input probability assignment for X E ctxJVc w ^ covar i ance Covjvec (X^ } 
Q T <g) In c subject to the average power constraint of the form 

^ E (tr(XXt) } = -L tr {Q T ® I Nc ) = tr (Q) < V (3) 

where Q is the tit x tit positive semidefinite matrix and V is the total transmit power over nx 
transmit antennas. Taking into account spatial fading correlation at both the transmitter and the 
receiver, we consider the channel matrix H is given by [17], [18] 

H = $i /2 # $T /2 (4) 

where $t £ c™ TXnT > and $r E C nRXnR > are the transmit and receive correlation 
matrices, respectively, and H ~ A/"„ Hi „ T (0, I nn ,I nT ) is a matrix with i.i.d., zero-mean, unit- 
variance, complex Gaussian entries. The entry H^, i — 1,2, ... , n R , j = 1,2, ... , ra T , of 
H is a complex propagation coefficient between the jth transmit antenna and the ith receive 
antenna with E {\Hij\ 2 } = 1. Note that H ~ J\fn K ,n T (0,$ R ,$ T ) [21]. With perfect CSI at the 
receiver, we have the transition pdf 

p(Y\X,H) = (irN o y n * N <etr{-±-(Y-HX)(Y-HX)^ (5) 

which completely characterizes a block-fading MIMO channel. 
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In what follows, we define the random matrix 9 G C mxm > as 

6 4 \ HH] i if n R ^ n T 6 
| H^H, otherwise 

which is a matrix quadratic form in complex Gaussian matrices, denoted by G ~ Q m ,n (In, $i, $2) 
[21], where m = min{nx,nR}, n = max {nx, hr}, and 

($ ieC ™$ 2 GC-) = | ($R ' $T)l ifn R <n T 

1 (<& T ,<I> R ) , otherwise. 

The pdf of 9 ~ Q m , n $ 2 ) is given by [22] 

Pe (e) = ^— det det($ 2 )- m det(er m F ( " ) (-$r 1 e,$ 2 1 ), 6 > 0, (8) 

T m (n) 

where f m (a) = 7]- m (™-i)/2 TJ™^ 1 r (a — i), K {«} > m — 1, is the complex multivariate gamma 
function, T (•) is the Euler gamma function, and p Fq a ^ (•) is the hypergeometric function of two 
Hermitian matrices, defined by [24, eq. (88)]. 

III. MIMO Exponent: Reliability-Rate Tradeoff 

This section is based on Gallager's random coding bound on the error probability of maximum- 
likelihood (ML) decoding for a channel with continuous inputs and outputs [3]. Notably, the 
bound determines the behavior of the error probability as a function of the rate and the codeword 
length. Hence, by determining Gallager's exponent, we can obtain significant insight into the 
reliability-rate tradeoff in communication over MIMO channels and the required codeword length 
to achieve a certain level of reliable communication. In particular, the diversity-multiplexing 
tradeoff of MIMO channels [12], [13], [26] is a special case of the reliability-rate tradeoff as 
the SNR goes to infinity. 

A. Random Coding Exponent 

Using the formulation developed in [3, ch. 7], we obtain the random coding bound on the 
error probability of ML decoding over block-fading MIMO channels as 3 

p e < (^—^ e -N b N c E r (p x (X),R,N c ) ^ 

3 When X = (Xij) is an m x n matrix of complex variables that do not depend functionally on each other, dX = 
UZi n"=i d&Xij dSXij . If X e C mxm is Hermitian, then dX = ]\™ =1 dX u f]^ dStX^ d^X tJ . 
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where r, 5 > and 



(10) 



4 



/2nN b al 

[ [ti(XX^-N c V] 2 Px (X)dX. (11) 
Jx 

The random coding exponent E x (p x (X ) , R, N c ) in © is given by 

E r ( Px (X) , R, N c ) = max (max E (p x (X) , p, r, N c ) - pR} (12) 

with 

E (p x {X),p,r,N c ) 

(13) 

The parameter r to be optimized may be viewed as a Lagrange multiplier corresponding to the 
input-power constraint [7]. 

1) Capacity-Achieving Input Distribution: As in [3]-[8], we choose the capacity-achieving 
distribution for p x (X) satisfying the power constraint ©, namely, 

Px (X) = 7r-" TjVc det {Q)~ Nc etr (-Q^XX^) (14) 

with tr (Q) < V. 4 Although this choice of the Gaussian input distribution is optimal only if the 
rate approaches the channel capacity, it makes the problem analytically tractable [3]. 

Proposition 1: Let E 0J ^- (Q, p, r, iV c ) be £"0 (p x (X) , p, r, N c ) in (fT3l for the Gaussian input 
distribution p x (X) of (TT~4T> . Then, we have 

E 0J t (Q, p, r, N c ) = rV (1 + p) + (1 + p) In det (I„ T - rQ) 

-^44- + * w ^+rr |- (15) 

Proof: See Appendix lAl □ 

4 In general, optimization of the input distribution px (-X') under the power constraint to maximize the error exponent (i.e., 
to minimize the upper bound) is a difficult task. 
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For the case of equal power allocation to each of transmit antennas, i.e., Q = —I nT (because 



the transmitter has no channel knowledge), (fT31) becomes 

E o,N (£Jnr, P, r, N c )=rV{l + p) + n T (1 + p) In 



N c ~~ | V'" R ' (n T - rP) (1 + p) 



-— lnE<(det ( J nR + - J, (16) 



where 7 = V/N is the average SNR at each receive antenna. Let us introduce a new variable 
/3 = n T — rV where /3 is restricted to the range < /3 < n T to have a meaningful result in (fl6l) . 
Then, we have 



E (p,(3,N c )± E ^^I nT ,p,r,N c 



1 



where 



= (1 + p) (n T - /3) + n T (1 + p) In (/3/n T ) -— In C (p, /?, AQ (17) 

V V ' 

C (p, p, N c ) = E I det (l m + - J® + p) ) "J. (18) 

With maximization over /3 G [0, n T ] and p G [0,1] to obtain the tightest bound, we have the 
random coding exponent for Gaussian codebooks and equal power allocation as follows: 5 



E r (R,N c )±E r ( Px (X),R,N c 



X~Af nT>Nc [0,^I nT ,I Nt . ) 



= ^\^ E ^^- PR )- (19) 

Proposition 2: Let (3* (p) be the value of (3 that maximizes E (p, /3, N c ) defined in (flTl) for 
each p G [0, 1]. Then, (3* (p) is the solution of 3E (p, /5, N c ) /d(3 — and is always in the range 

< /? < n T . 

Proof: See Appendix |BJ □ 
It can be shown using (l64l) and d67t in Appendix |B] that as 7 — > 00 or 7 — ► 0, the optimal 
value of [3 does not depend on N c , that is, 

lim (3* (p) = n T ' n W_ an( j ^ m ^* ^ __ n ^ 

7^00 1 + P t - * 

5 The random coding bound can be improved by expurgating "bad" codewords from the code ensemble at low rates (see, e.g, 
[3]). More details for the expurgated exponent of block-fading MIMO channels can be found in [27]. 
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According to Proposition |2] and using the general relation dE r (R, N c ) jdR = — p, the maximiza- 
tion of the exponent in (fl9t over (3 E [0,71t] and p E [0, 1] can be performed by the following 
parametric equations: 

E T (R,N c )=E (p,{3* (p),N c )-pR 
'dE (p,(3,N c 



R 



dp 



(20) 
(21) 



0=0* (p) 



with 



dE (p,(3,N c 



dp 



(n T - 0) + n T In ((3/n T )-—C^ (p, (3, N c ) 



dC (p,P,N c 
dp 



(22) 



where 

dp 



-N c p 



PI 



0(i+pY 



tr^e(4n^)-U-indet(4n^) 



(23) 



2) Key Quantities: The values of R in d2TT) at p = 1 and p = are the critical rate R CI and 
the ergodic capacity (C) of the channel, respectively [3]-[6]. From dE (p, (3, N c ) /d(3 in (l67l) . 
we see that /3* (0) = and hence, the ergodic capacity can be written as 

~dE (p,(3,N c 



(C) 



dp 

E jlndet (l m + 



p=0, 0=riT 



— 



(24) 



(25) 



We remark that the capacity expression d23t obtained from the exponent is independent of the 
channel coherence time N c and is in agreement with the previous result [14]— [16]. Also, the 
quantity E is defined as the value of the exponent E r (R, N c ) at R = 0, referred to as the 
exponential error-bound parameter [4], [5], and is given by E (1, (3* (1) , N c ). This quantity is 
equal to the value of R at which the exponent becomes zero by setting p = 1 and (3 = (3* (1). 
If setting r = or equivalently (3 = n T (i.e., without the constraint on the minimum energy of 
the codewords) in (fT3l . E becomes equal to the cutoff rate R of the channel 

(26) 



R = E {l,n T ,N l 
1 



InE <^ det [I m + 



2nT 







-Nc 



(27) 



This is an important parameter, as it determines both the magnitude of the zero-rate exponent 
and the rate regime in which the error probability can be made arbitrarily small by increasing 
the codeword length. 
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3) Effect of Channel Coherence — Channel-Incurable Effect: Using Jensen's inequality, it is 
easy to show 

ln£ (p,/3,7V c -l) (28) 



^ln£ (p,AiV c )>^i 



1 



yielding 



E (p,(3,N c )<E (p,(3,N c -l). (29) 

Therefore, for fixed R, the random coding exponent decreases with iV c , while the channel capacity 
is independent of N c . This reliability reduction is due to the fact that the increase in N c results in 
a decrease in the number of independent channel realizations across the code and hence, reduces 
the effectiveness of channel coding to mitigate unfavorable fading. We call this effect of the 
channel coherence time on communication reliability "a channel-incurable effect". In particular, 
since lim^^oo In C (p, @, N c ) = 0, we have 

lim E Q {p,P,N e )=K.{p,P) (30) 

A< c -+oo 

leading to limjv^oo P* (p) — and limjv^oo E T (R, N c ) = 0. Therefore, if iV c — > oo, it is 
impossible to transmit information at any positive rate with arbitrary reliability even with the use 
of multiple antennas. In fact, n T must also increase without limit so that the so-called space-time 
autocoding effect takes place, which makes arbitrarily reliable communications possible [11]. 

B. Evaluation of E (p, (3, N c ), dE (p, (3, N c ) /d(3, and dE (p, (3, N c ) /dp 

To calculate the random coding exponent, the quantities E (p, (3, N c ), dE (p, (3, N c ) /d(3, and 
8E (p, (3, N c ) j dp need to be determined. We now evaluate them in the following theorem which 
will be proven in the next section. 

Theorem 1: Let H ~ A/" nRi „ T (0,* R ,* T ) or 6 ~ Q m>n (7 n ,*i,* 2 )- Then, 

1) E (p,f3,N c ) is given by 



In [K-}det 



if N c p E {l,2,...,m} 



E (p,f3,N c ) = { 



-^ln \%{p,N e ) det 



G(n-m) (^2) 

r( P ,(3) 



otherwise. 



(31) 
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If $ T = I riT and $ R = J„ R (i.i.d. MIMO channel), then E (p, P, N c ) reduces to 



E (p, P, N c ) =JC(j>,/3)-— lnf^ 1 det T iid (p, /3) ) . 



2) d£ 0, P, N e ) /dp is given by 



dE (p,l3,N c ) W/JW m T A I 



dp 



1C W) (P,P)~ 



- — tr 



•(n-m) ($2 

T(p,0) 



-1 r 





rw (p, /3) 



If $ T = J nT and $ R = I nR , then 

gte^j = K<» („, P) - ± tr {T i> (p, T«« („,/?)} 
3) dE (p, p, N c ) I dp is given by 



dE Q (p,/3,N c ) 
dp 



(p, p) 



m— 1 



E 



iV c (l+p) 



tr 



(32) 



(33) 



(34) 



v y i 

pl N c p-m + Xi (*i) - J 



^N c p-k N c 



'{n-m) ($2j 

T(p,/3) 



-I r 





TW (p,/3). 
N c p^l,2,...,m-1. (35) 



If $ T = I riT and $ R = J na , then 

^y^ c) = (p, /?) - ^ tr {T- (p, /3) T« (p, /?) } . (36) 

The quantities K cm , K iid , T A , T B (p, N c ), and the matrices G ( .) (•), 2 (p,P), T (p,/3), T (/3) (p,/3), 
TW (p, T iid (p, /3), T<J (p, P), and (p, /?) are given in Table ffl 

Corollary 1 (Ergodic Capacity): If H ^ A/"„ Ri „ T (0,$r,$ t ), then the ergodic capacity (C) is 
given by 



(C) = tr 



G(n-m) ($2) 

v _ T(0,n T ) 
with A G R mxn given by 



m 



e(*i) x*(*i)-i 

!) + E E — T 



(37) 



Ai,i • • • Ai (e (* 2 ) 
A^*!),! • • • A e( * l)ie (* 2 ). 



(38) 
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where the (i,j)th entry A p>t}jij of A p , q £ RXp(*i)xx 9 (* 2 ) 5 p = i ? . . . , e g = 1, . . . , e ($ 2 ), is 
A P ,™ = \ P) A (9) ($ 2 ), m - * + l) . (39) 

Proo// Note that 



tr 



G(„_ m ) ($2) 
T(0,n T ) 



-1 







TW (0,n T ) 



-iNLtr 



+ nx tr 



G(n-m) ($2) 


-1 


"0" 


) 


. T(0,nr) . 




A 




G(n-m) ($2) 


-1 









_ T(0,n T ) _ 




T (/3) (0,n T ). 


} 



(40) 



-r A 



The proof follows immediately from (|24b . Theorem [T] 3 with p = and /? = tit, and d40l . □ 
Note that the expression (071) for the ergodic capacity (C) is sufficiently general and applicable 
to arbitrary structures of correlation matrices $ T and $ R , and hence, embraces all the previously 
known results for i.i.d. channels ($x = In T , $r = I nR ) [16], [21], one-sided correlated channels 
($1 = I m [19] or $ 2 = I n [20]), and doubly correlated channels [22] (where all the eigenvalues 
of $ T and $ R are assumed to be distinct) as special cases of (l37t . 

Corollary 2 (Cutoff Rate): If H ~ A/" nRjnT (0,$ R ,$ T ), then the cutoff rate R is given by 



-^m I K~} det 



#0 = < 



G( m _AT c ) ($1) 

S(l,nr) 



f ln (^)-A ln ( T B(l,iVc)det 



G(n-m) ($2) 

T(l,n T ) 



if JV C G {1,2,..., m} 



, otherwise. 



(41) 



In particular, if $ T = 7„ T and $ R = J nR , then we have 

R = ~ ln(^ det T iid (1, tit)) • (42) 
Proof: It follows immediately from (1261) and Theorem [T]l with p = 1 and (3 = n^. □ 

C. Coding Requirement 

As in [6], we can approximate the required codeword length to achieve a prescribed error 
probability P e at a rate i? by solving for N b in the following equation: 



P e =(2e rS /^)e 



2 „-N b N c Er(R,N c ) 



(43) 
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Using (flOb . it is easy to see that the factor (2e r5 /£) 2 in (l43l is minimized over 5 > 0, for large 
N^, by choosing § = 1/r [3]. This yields 

min (2e rS /g) 2 « 8vreV?r 2 A^ b for large JV b . (44) 

Also, from (fTTT) and [23, Lemma 5], we have 

a | = N c V 2 /n T . (45) 

Combining (l44l) and (|431) together with the fact that (3 = — rV, d43l can be written as 

P e = (8vr/n T ) {n T - p* (p)f N h N c e -N b N c E t (R,N c )+2 _ (46) 

After solving for Ar b in d46b . we take L = N c ■ \N b ~\ as our estimate of the codeword length (in 
symbol) required to achieve P e at the rate R, where [•] denotes the smallest integer larger than 
or equal an enclosed quantity. 

IV. Proof of the Main Theorem 

In this section, we provide proofs of the main results stated in Theorem [T] The methodology 
recently developed in [22] and [23] for dealing with random matrices paves a way to prove the 
theorem. 

A. Proof of Theorem [7] 1 

Using the same steps leading to [22, Theorem 1], we get 

C (p, 0, N c ) = [ det (J m + rfi)-"* pe (6) dQ 
</e=et>o 

l: dct ($,) '" ' "' 



, , n A r m ( )n(^( e )- A ,( )) 2 

T m in) T m [m) A(e) 

x ^ (Kp;D, -ntd F {n) (D, dX (6) (47) 

where rj = ^r^y and D = diag (A x (G), A 2 (0), . . . , X m (O)). 6 Successively applying the generic 
determinantal formula for hypergeometric functions of matrix arguments [23, Lemma 4] and the 

6 For A = A* £ C pxp > 0, we denote J X(A) dX {A) = fo*- l(A) J^A)^ ■ ■ ■ J^ {A) dAj {A)dX 2 (A) ■ ■ ■ d\ v {A). If the 
integrand is symmetric in \ 1 (A), X 2 (A), . . . ,X p (A), then 



dA (A) = i / • • • / dAi (A)dA 2 (A) ■ ■ ■ dX p (A). 
X(A) P- Jo Jo 

p-fold 
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generalized Cauchy-Binet formula [22, Lemma 2], the integral in (l47t can be evaluated, after 
some algebra, as 



C (p,P,N c )= V - T -T B (p,N c ) det 



G(n-m) ($2) 

T(p,/3) 



iV c p^l,2,...,m-l. 



(48) 



Substituting d48t into (fTTt gives the second case of (I3TT) . It should be noted that the formula 
in the second case of (I3H has singular points at iV c p = 1, 2, . . . , m — 1 for each p 6 (0, 1]. These 
singularities stem from the quantity T B (p,N c ), which can be alleviated using the following 
analysis. 

Suppose that N c p is a positive integer. Then, using [23, Lemma 1], we have 
£ (p,/3,i\Q=Ee{E s {etr (-776^)}} 

= E s {det (I mn + rjSS^ ® $ 2 ) _1 } 

= E s {det (7 mn + ^ t $ 1 ^®$ 2 )" 1 } (49) 

where 5 ~ M m ,N cP (0, I mi lN cP ) is a complex Gaussian matrix statistically independent of G, and 
the last equality follows from the fact that SS^ 1 and S^iS have the same nonzero eigenvalues. 

If N c p E {1,2, . . . ,m}, then Z = S^iS ~ Q Nc p, m {I m ,I NoP ,Qi)- Hence, using [23, Theo- 
rem 9], (l49l for the case of N c p G {1,2,..., m} can be written as 



£ (p, /?, iV c ) = E x[z) I 11 det{ J„ + V \ k (Z) $ 2 

lfe=i 



-1 



K~t / H det \ln + Vh(Z)*2 

JX{Z) k=1 



x det 



G(m-AT cP ) ($1} 



J det (Ar 1 (Z)) d\ (Z) (50) 

/ l<i,j<N c p x 3 ' 



where a = 
is given by 

Now, applying [22, Lemma 2] to d50l yields 



and the (i, j)th entry H Mi of H fe e l^^' 1 * 1 ', k = 1, 2, . . . , g 



(51) 



£ (p, /3, N c ) = K'l det 



E(p,/3) 



(52) 
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where the (z, j)th entry E k ^ (p, (3) of the fcth constituent matrix H fc (p, f3) is given by 

/■oo 

E ktij (p, (3) = / det (J n + t^)" 1 z^-V^W** (53) 
Using the characteristic coefficients [23, Definition 6], d53l can be written as 

f?(*2)Xp(*2) ^ 

E k>ij (p,f3)=J2 E ^( $2 ) / (l + ^^aJzJ^^-V^w^tfa (54) 

p=l g=l ^° 

where X p>q ($ 2 ) is the (p, g)th characteristic coefficient of $ 2 . Finally, substituting (l52l into (fT7l 
gives the first case of d3"TT ) and hence, we complete the proof of the first part. 

B. Proofs of TheoremU\2 andU\3 

The second and third parts can be obtained by differentiating E (p, ft, N c ) in Theorem [T]l 
with respect to (3 and p, respectively, with the help of the logarithmic derivative of a determinant 
[28, Theorem 9.4] (or more generally [22, Lemma 1]). 

V. Numerical Results and Discussion 

In this section, we provide some numerical results to illustrate the reliability-rate tradeoff 
in block- fading MIMO channels. For spatial fading correlation, we consider an exponential 
correlation model with $t = (Ct J ') an d $R = (Cr _J ')> Ct> Cr [0, 1), in all examples. 

To ascertain the effect of the channel coherence on the error exponent, Figs. |2] and |3j respec- 
tively, show the random coding exponent E T (R, N c ) as a function of a rate R for i.i.d. (Ct = 0, 
Cr = 0) and exponentially correlated (£t = 0.5, Cr = 0.7) MIMO channels at 7 = 15 dB, where 
rvj = tir = 3 and N c ranges from 1 to 10. We can see from the figures that the exponent at a rate 
R below the ergodic capacity decreases with iV c , while the ergodic capacity remains constant for 
all N c (i.e., 8.48 and 7.19 nats/symbol for Figs. El and [3 respectively). For example, the error 
exponents at rates R < R CT for iV c = 10 are reduced by roughly 3.46 and 2.86 for i.i.d. and 
exponentially correlated cases, respectively, compared with those for N c = 1. This reduction in 
the exponent, consequently, requires using a longer code to achieve the same error probability. 
Hence, we see that unlike the capacity (with perfect receive CSI), the channel coherence time 
plays a fundamental role in the error exponent or the reliability of communications. 

Fig. E| demonstrates the effect of spatial fading correlation on the random coding exponent, 
where £t = ( R = (, 7 = 15 dB, n T = n R = 3, N c = 5, and ( ranges from (i.i.d.) to 0.9. As 
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seen from the figure, there exists a remarkable reduction in the exponent at the same rate due to 
correlation, especially for ( > 0.5. The amount of reduction in the exponent at rates R < R cr , 
relative to the i.i.d. MIMO exponent, ranges from 0.07 for £ = 0.2 to 2.17 for £ = 0.9, indicating 
that a longer code is required to achieve the same level of reliable communications. Equivalently, 
a decrease in the information rate is required for more correlated channels to achieve the same 
value of the exponent. For example, the exponent at a rate 3 nats/symbol are 1.94 and 1.53 for 
the i.i.d. and correlated (C T = Cr = 0.5) channels, respectively. This implies that 27% increase in 
the codeword length, due to spatial fading correlation, is required to achieve a rate 3 nats/symbol 
with the same communication reliability. 

To get more insight into the influences of the number of antennas, channel coherence time, and 
fading correlation on a coding requirement for MIMO channels, the codeword length required to 
achieve P e < 10~ 6 at a rate 8.0 bits/symbol (5.55 nats/symbol) are investigated in Tables HTl - IIVI 
The codeword lengths in the tables are calculated in such a manner as described in Section IIII-CI 
Table [H] serves to demonstrate the effect of increasing the number of antennas on the coding 
requirement, in which the required codeword length L is shown for i.i.d. MIMO channels with 
N c = 5. Note that it is impossible to reliably communicate at a rate 8.0 bits/symbol below the 
SNR 7 of 14.55 dB, 9.68 dB, and 6.79 dB for nx = = 2, 3, and 4, respectively, since these 
SNR's are required to attain the ergodic capacity (C) of 8.0 bits/symbol in each of the cases. As 
seen from the table, with increasing the number of antennas at both transmit and receive sides, 
the required codeword lengths are remarkably reduced. This is due to the advantages of the use 
of multiple antennas, e.g., spatial multiplexing and diversity gains [21]. For example, at 7 = 16 
dB, increasing the number of antennas at both sides from 2 to 3 and 4 reduces the corresponding 
codeword length to almost 2.8% and 0.9% of the amount required for two transmit and receive 
antennas, respectively, which is a tremendous reduction in the codeword length. 

Table [in] shows the required codeword length L for i.i.d. and exponentially correlated (£t = 
0.5, Cr = 0.7) MIMO channels with n T = n R = 3 at 7 = 15 dB when iV c varies from 1 to 10. It 
is clear from Table ITTT1 that for each value of N c , the codeword lengths for correlated channels are 
much longer than those for i.i.d. channels. For example, the increase in the required codeword 
length, due to exponential correlation (£ T = 0.5, Cr = 0.7), ranges from 194% for iV c = 1 to 
138% for iV c = 10, which is a significant increase in required codeword length. Also, when 
going N c from 1 to 10, there is a considerable increase in the required codeword length, relative 
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to that for the single-symbol coherence time, which ranges from 33% to 344% for the i.i.d. case 
and from 28% to 258% for the correlated case, respectively. 

Table |IV| demonstrates the effect of correlation on the required code length L, where n T = 
% = 3, Cr = Cr — C> N c = 5, and 7 = 15 dB. The table contains the corresponding codeword 
lengths for £ from to 0.9. As seen from the table, the required codeword length for the 
case of exponential correlation £ = 0.7 is equal to 4.5 times as long as for the i.i.d. channel 
(C = 0). Particularly, when £ > 0.5, there exists a large amount of increase in required codeword 
length due to a stronger correlation. Also, since the ergodic capacity is 7.36 bits/symbol for 
Ct = Cr = 0.9 at 7 = 15 dB, it is impossible to achieve reliable communications at a rate 8.0 
bits/symbol (regardless of the codeword length), when Ct = Cr = 0.9. 

Finally, Fig. |5] shows the cutoff rate Rq in nats/symbol as a function of a correlation coefficient 
C for exponentially correlated MIMO channels with £ T = £ R = £ at 7 = 15 dB, where n R = 
fiR = 3 and N c varies from 1 to 10. We see that the cutoff rate Ro decreases with N c for all 
C £ [0, 1). While (C) remains constant, R monotonically decreases with N c , going to as 
N c — > 00 (see (l29l and (l30l with p — 1 and (3 = n T ). Hence, these two measures diverge as iV c 
increases and eventually lirn/v c -*. 00 7^ = 00, which coincides with the divergent behavior of the 
capacity and cutoff rate of a channel with block memory [9]. This observation reveals that R 
is more pertinent than (C) as a figure of merit that reflects the quality of block-fading channels. 

VI. Conclusions 

In this paper, we derived Gallager's random coding error exponent to investigate a fundamen- 
tal tradeoff between the communication reliability and information rate in spatially correlated 
MIMO channels. We considered a block-fading channel with perfect receive CSI and Gaussian 
codebooks. The required codeword lengths for a prescribed error probability were calculated from 
the random coding bound to aid in the assessment of the coding requirement on such MIMO 
channels, taking into account the effects of the number of antennas, the channel coherence time, 
and the amount of spatial fading correlation. In addition, we obtained the general formulae for 
the ergodic capacity and cutoff rate, which encompass all the previous capacity results as special 
cases of our expressions. In parallel to the capacity-cutoff rate divergence in a block- memory 
channel, we observed the channel-incurable effect: the monotonically decreasing property of the 
MIMO exponent (i.e., communication reliability) with the channel coherence time. 
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Appendix 

A. Proof of Proposition \J\ 

Lemma 1: Let S ~ N m , n (M,E, J n ) and A e C mxm > be Hermitian. Then, we have 

E {etr (-ASS r ) } = det (J m + HA)~ n etr {- (A 1 + E) ' MM 1 } . (55) 
Proo/; Note that 

E {etr (-ASS 1 ) } = (1<<r ^ I etr { - ASS 1 - E 1 [S - M) (S - M) 1 } (IS. (56) 

By writing the trace of the quadratic form in the exponent of (1561) as 

tr {ASS* + E- 1 (5 - M) (5 - M) f } 

= tr{(A I E ') [S- (/„, + T.A ) 'M] [S- (/,„ +EA) 'M] f - (A 1 +E) 'MM 1 }, 

(57) 

we get 

E {etr (-AS&) } = det ^} - etr {- (A ' + E) ' MM 1 } 

x y etr {- (A + E ' ) [5 - (/„, + HA) 1 M] [5 - (/„, + T.A) 1 M] '} r/S (58) 

v ' 

= n mn dct^+E- 1 )"™ 

from which (l55l) follows readily. □ 
Proof of Proposition UJ Using Lemma [H we have 



/ (X) e r H xx ')- N ^] p (Y\ X, H) l/i - l+p) dX 
Jx 



= e"^ (TriVo)-^^ 1 ^ det (J„ T - rQ)^ det (l nR + g ^ 

\ ^ V (-L + PJ 

x etr f- / - (j BH + H^-rI n ^W\ ^ ^ 

\ N (l + p) \ nR N (l+p) J 1 

Substituting (l59t into (fl3l) and integrating over V, we have 

J p x {X)e r H xx 'y N <%(Y\X,H) l/ ^Uxy P dY 

= e -^d + p) det (JnT _ rQ) -iVcd + p) det fj nn + H (Q-^-rl^- 1 # ] (fiQ) 

Finally, substituting d60l) into (fT3t yields the result (fl5l) . 



-JV c /> 
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B. Proof of Proposition 

We provide a sketch of the proof of Proposition |2] using a similar approach in [4] and [6]. 
For notational simplicity, let us denote Clp^ = (31 m + 76/ (1 + p). Then, E (p, (3, N c ) in (flTt 
can be rewritten as 

E (p, (3, N c ) = K (p, /3) - mp In {(3) - ±- In L x (p, /3, iV c ) (61) 

where £ : (p, (3, N c ) = E{det (n^) - ^}. Since K (p, /3)-mp In (J3f) is concave in /3, E (p, /3, JV C ) 
is a concave function of /3 if InC^ 1 (p, /?, iV c ) is concave in (3 for all p G [0, 1]. Noting that 

92 A - (p, ft JVJ { (M4gm 2 - £, (p, ft JV„) 92£ ' A * 



and £ 1 (p, (3, N c ) > 0, it is sufficient to show that 



(62) 



a^MNjy p^mn^ (63) 



It is easy to show that 

^|W =d et(n p ,)t r (n-a and £^ = -tr(fl 3 ) 

and hence, 

^ iP d ^ Nc) = E{-N c p det (fl^)-^ tr (n-y } (64) 
d^P^Q = E det (n^)-^ [iv c p tr 2 (fi-J) + tr (O; 2 )] } . (65) 
Let us now define the random variables 

X 2 = det (Q p ,p)- N ' p and Y 2 = (7V c p) 2 det (0^)"^" tr 2 . 

From Schwartz's inequality, we have 

2 



M?JMlM\ = E 2 {XY}<E{X 2 }-E{Y 2 } 



< E {X 2 } ■ E {Y 2 + N c pX 2 tr (fi~ 2 ) } 
dp 2 

From (R^oT). we see that E (p,(3,N c ) is a concave function of (3 for all p e [0, 1]. Hence, the 
maximum over (3 occurs at (3* (p) for which [8E (p, (3, N c ) /d(3] \ f3=f5 ,^ p - ) = and it is sufficient 
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to show that [dE (p, 0, N c ) /OP] \ p=0 >0 and [dE (p, p, N c ) /dp] ^ < for all p G [0, 1] 
in order to prove < P* (p) < n T . Since 



dE (p, P, N c ) _ (1 + p) (n T -P) mp 1 . ^ dC, (p, p, N c ) 

dp P P N c (IK >j dp 

V v ' 



it is clear that lim^o dE (p, /3, iV c ) — > oo. Also, 



dEo(p,P,N c ) 
dp 



/3=n T 



m ^ 1 AM 



dp 



(3=n T 



Since £ x (p, /3, JV C ) > and 



£ x (p,n T ,iV c ) =E{det (O p , nT )-^} = E jdet (fi p , nT ) 
>^E{det(O p , nT )-^ tr(fl-y} 



Ut 1 

mp iV c 



(p,P,N c 



dp 



0=n T 



it follows that 



■^-^r 1 (p,n T ,JV c ) 
n T iV c 



dp 



< 0. 



f3=1lT 



Thus, [9^0 (p, iV c ) /<9/3] < and we complete the proof of the proposition. 



(67) 



(68) 



(69) 



(70) 
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In Theorem 

G w (*) = [^tj A 2 ■•■ i #) ]eK" xp , for # p x p Hermitian, v<p 

G M (*)= [Ax a 2 ••• A em ] ew xp 

where ' 

A k = {Ak,ij) ei"^'*', * = i,2,..., e (*) 

element: A,,, = (-1)'^ {i - j + l),_ a X^ +j (») 

A = 6 K" XXfe W , fc = 1, 2, . . . , q (») 

(i, j)^ element: 34 Mj = M)^ (< - j + A*"/ (*). 

^ ( a )n = h (a + 1) ■ • • (a + n — 1), (a) = 1 is the Pochhammer symbol. 



/c COT = dot($ 1 ) jVtf, dct{G (m) ($i)} n (*- 


-1)!, 


N e pe{l,2,...,m} 


Km = II (» - ( fc - 1)! 

fc=i 


e(*i) r 

T A = ±m(m + E »(*i)x<(»i) 

i=i L 






T B (p,iV c ) =det(* 2 )- m det{G (m) 


det {G(„) 





B(p,/3) = [s 1 (p,/9) 3 2 (p,/3) ■■■ B e(#l) 0>,/3)J el"' pxra 
where 2 ' 

E fc (p, 0) = (H fc ,^ (p, 0)) e R^x*>=(*i>, fc = 1, 2, . . . ,g («!), AT p e {1, 2, . . . , m} 
(i, j)— element: 

Sk,« (P./3) = E^L* 2) E*=i* 2) (*a) A (p> ($2), A (fc) («!), -g + 1 

2) A'p,, ($2) is the (p, g)fh characteristic coefficient of $2 (see for details [23, Definition 6]). 
The Q KtV (a, b, p) is defined as the integral 

00 

Q K . V (a, b, p) = jf (1 + ax) 1 *- 1 In"- 1 (1 + ax) x K ' 1 e - x/b dx, a, b > 0, «, 1/ G N b , /U G C 


6" (k — 1)! 2J0 (k, — M + 1; — afc) , if v — 1 



a- (./ - 1)! e 1 ^ 6 ' E 1 [(-I)"-*" 1 (V) (abr +k Gft&^fr^,^)] , 



otherwise 



where p F q (01, 02, • • • , dp; 61, 62, • • • , b q ; z) is the generalized hypergeometric function of scalar argument 
[29, eq. (9.14.1)] and G™ 9 ™ (■) is the Meijer G-function [29, eq. (9.301)]. The detailed derivation of this 
integral identity can be found in [22, Appendix A]. 
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TABLE I 

(Continued.) SOME QUANTITIES AND MATRICES INVOLVED IN THEOREM^ 



In Theorem Q] 



T(P,£) 



T M (p,/3) 



T „(«!),! (p,/3) 



Ti lfi (# a ) (p, /3) 



• e(*i),e(#2) 



where 



T M (p, 13) = (T p , q , I3 (p, P)) G M^(*i)xx <f (* 2 ) ; p = l, 2 , .... e («!),« = 1, 2, .... e ($ 2 ) 

element: T p , q ,ij (p, /3) = Qi+j-i,i ( 0( i +p) ^(p) {&i),\ q ) (#2), -iV c p + m - i + 1J . 



T<« ( Pj /3) 4^T(p,/9) = 



T M(* 2 )(P./3) 



."^e(*i),l (Pj^) '" ^e(*i),e(*2) (P'@\ 

*l)xx g (*2) 



where 

T(«(p,/3) = (T^.(p,/3)) g: 
(i, j)— element: 

T p 9 ? } ,v (P" 9 ) = ^(i+p) (^Vcp-m + i)A (p) (^y^jA^) (*i),A <9> (* 2 ), -iV c/9 + m • 



TW (p,/?)^|T(p,/?) = 



tW(p,/3) 



•l)xx«(«2) 



i,e(*2) 



(p,/3) 



where 

T^(p,/3) = (TW >i .( P) /3)) G 
(i, j)— element: 



T Ud (p,p) = {T m , ij (p,/3))GR mxm 

(z,j)- element: Tad.ij (p,/3) = t/n-m+i+j-i.i ( g (1 7 +p) » 1. ~N c p + 1 



-Mid 


(p,0)±-feTm (P,I3) = 
element: 


iid , 

(p,/3) 


(p,/3)) G R mxm 


, 1, -JVep) 




_ /3 2 (l+p) yn-m+t+3,1 yp( 1+p ) 


-y-(p) 
Mid 


element: 


iid f ij 

(p,/3) 


(p,(3)) G R mxm 

= t^ t ^(p^) - N c g n . m . 


H+i-l,2 (/J(l+ rt .l. ^cP+l) 
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TABLE II 

Required codeword length L as a function of SNR 7 for i.i.d. MIMO channels (<fr T = I nT , $r = I nR ) at a 

RATE 8.0 BITS/SYMBOL WITH P e < 10~ 6 FOR DIFFERENT NUMBERS OF ANTENNAS AND N c = 5 



SNR (dB) 


Codeword length L 


n T = n R = 


2 


n T = TiR. = 3 


tit = % = 4 


8 








510 


10 






10865 


75 


12 






210 


30 


14 






65 


15 


16 


1070 




30 


10 


18 


205 




20 


5 


20 


90 




15 


5 



Note: The ergodic capacity (C) of 8.0 bits/symbol is attained at 
7 = 14.55 dB for n T = n R = 2; 7 = 9.68 dB for n T = n R = 3; 
and 7 = 6.79 dB for tit = n R = 4, respectively. 



TABLE III 

Required codeword length L as a function of channel coherence time N c for i.i.d. and exponentially 

CORRELATED MIMO CHANNELS AT A RATE 8.0 BITS/SYMBOL FOR P e < 10~ 6 , n T = n R = 3, AND 7 = 15 dB 



Coherence time N c 


Codeword length L 


i.i.d. 


Ct = 0.5, Cr = 0.7 


1 


18 


53 


2 


24 


68 


3 


30 


84 


4 


36 


100 


5 


45 


115 


6 


48 


126 


7 


56 


140 


8 


64 


160 


9 


72 


171 


10 


80 


190 
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TABLE IV 

Required codeword length L as a function of correlation coefficient ( for exponentially correlated 
MIMO CHANNELS WITH Ct = Cr = C AT A RATE 8-0 BITS/SYMBOL FOR P e < 1CT 6 , n T = n R =3, N c = 5, AND 7 = 15 dB 



Correlation coefficient £ 


Codeword length L 


0.0 


45 


0.1 


45 


0.2 


45 


0.3 


50 


0.4 


60 


0.5 


75 


0.6 


105 


0.7 


200 


0.8 


1015 


0.9 





Note: For (t = Cr = 0-9, the ergodic capacity (C) is 
7.36 bits/symbol at 7 = 15 dB. 
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Fig. 1. A wireless communication link with tit transmit and riR receive antennas to communicate at a rate R over JVb 
independent iV c -symbol coherence intervals. 
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R (nats/ symbol) 



Fig. 2. Random coding exponent for i.i.d. MIMO channels (£ T = 0, Cr = 0) when N = 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10. 
nr = iir = 3 and 7 = 15 dB. 
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Fig. 3. Random coding exponent for exponentially correlated MIMO channels with £/r = 0.5 and £r = 0.7 when N c — 1, 

2, 3, 4, 5, 6, 7, 8, 9, and 10. n T = n R = 3 and 7 = 15 dB. 
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Fig. 4. Random coding exponent for exponentially correlated MIMO channels when £t = Cr = (i.i.d.), 0.2, 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8, and 0.9. n T = n R = 3 and 7 = 15 dB. 
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Correlation coefficient, C, 



Fig. 5. Cutoff rate Ro in nats/symbol as a function of a correlation coefficient £ for for exponentially correlated MIMO 
channels with ( T = Cr = C when iV c = 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10. n T = n R = 3 and 7 = 15 dB. 
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